Shot Noise in Ballistic Quantum Dots with a Mixed Classical Phase Space 



H.-S. Sim and H. Schomerus 
Max- Planck- Institut fur Physik komplexer Systeme, Nothnitzer Str. 38, 01187 Dresden, Germany 

(Dated: February 1, 2008) 

We investigate shot noise for quantum dots whose classical phase space consists of both regular 
and chaotic regions. The noise is systematically suppressed below the universal value of fully chaotic 
systems, by an amount which varies with the positions of the leads. We analyze the dynamical origin 
of this effect by a novel way to incorporate diffractive impurity scattering. The dependence of the 
shot noise on the scattering rate shows that the suppression arises due to the deterministic nature 
of transport through regular regions and along short chaotic trajectories. Shot noise can be used to 
probe phase-space structures of quantum dots with generic classical dynamics. 
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The phenomenon of shot noise, the time-dependent 
fluctuations in electrical currents caused by the discrete- 
ness of the electron charge e, has been extensively stud- 
ied in recent years in mesoscopic systems (for a review 
see ) . The focus of interest was on classically chaotic 
ballistic quantum dots (electron billiards) and on dif- 
fusive quantum wires, which can be investigated the- 
oretically by applying random matrix theory J2|, |^, ^] 
as well as semiclassical methods |$], ^J, [7j. For incoher- 
ent transport through a chaotic quantum dot the shot 
noise would assume the Poissonian value Pq = 2cGqV, 
where Go = Ne 2 /(2h) is the serial conductance of the 
two quantum point contacts which connect the dot to 
electronic reservoirs (maintained at a voltage difference 
V) by N channels. For low temperatures correlations of 
electrons due to Fermi statistics suppress the noise P by 
a factor of T = Pj Pq relative to this value of uncorre- 
cted electrons. In chaotic quantum dots the suppression 
factor T = J- c \^ = 1/4 is universal |3J, i.e., indepen- 
dent of the details of the system, which also has been 
confirmed by an experiment [g. The origin of the low- 
temperature noise is the probabilistic nature of quantum 
transport, arising from attempts to transmit charge car- 
riers between electronic reservoirs with a finite success 
probability <= [0,1]. Very recently, non- universal correc- 
tions to the shot-noise due to residual signatures of clas- 
sically deterministic scattering |H] have been discussed 
by Agam, Aleiner, and Larkin pX] . They predicted that 
shot noise in a chaotic dot can be further reduced be- 
low T c h under the condition that electrons pass the dot 
without sufficient diffraction. This has been verified in a 
recent experiment fll|| . 

In this paper we address the shot noise of generic bal- 
listic quantum dots. They are not classically chaotic Q 
but possess a mixed phase space, where regular islands 
are separated from chaotic seas by impenetrable dynami- 
cal barriers. Signatures of the mixed phase space in quan- 
tum transport have been found in the conductance, which 
exhibits fractal fluctuations (ll| and isolated resonances 
[pj. Shot noise can be seen as the second cumulant of 
charge counting statistics, with the conductance being 




first cumulant. As we will demonstrate, shot noise carries 
valuable dynamical information which can be extracted 
systematically. For generic quantum dots the suppres- 
sion factor T is found to be systematically reduced be- 
low the universal value T&i for fully chaotic systems, not 
only due to the deterministic nature of transport along 
short chaotic trajectories (this mechanism of Ref. JjC[ ] 
will be confirmed), but also because quantum diffraction 
is strongly reduced for transport through regular regions, 
as well. 

Our investigation of the dynamical origin of the addi- 
tional shot-noise reduction is based on a novel procedure 
to analyze shot noise (which we obtain from a numerical 
simulation) with help of the Poisson kernel pH 
a statistical ensemble of random matrix theory 
the Poisson kernel one averages the scattering matrix 
over an energy range E av , in this way eliminating the 
system-specific details of the dynamics with time scales 
longer than i av = h/E avl and replaces these by random 
dynamics of the same universality class as elastic diffrac- 
tive impurity scattering (equivalently, fully chaotic dy- 
namics). (Ref. [[Tot also employed diffractive impurity 
scattering; however, there it was used to mimic the non- 
deterministic processes of quantum chaos within a theory 
that cannot account for them, while in our case all quan- 
tum effects are fully included from the beginning.) The 
effective mean free scattering time £ av can be tuned by 
changing the energy-averaging window i? av . The depen- 
dence of the suppression factor T on t av then allows to 
analyze the properties of trajectories with classical dwell 
times tdweii — *av (up to a possible factor of order unity) , 
over the large range of dwell times that is typically in- 
volved in the classical transport through generic quantum 
dots (in contrast, chaotic dots are characterized only by 
a single time scale, the mean dwell time). The analysis 
via the Poisson kernel is supplemented by a semiclassical 
estimation of the shot noise, obtained from classically de- 
terministic motion course-grained over Planck cells while 
preserving Fermi statistics, in the spirit of Refs. [^,0- It 
follows that shot noise is suppressed most strongly if the 
leads are well-coupled to classical regular regions (with 
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FIG. 1: (a) Schematic diagram of the annular billiard be- 
tween two circular hard walls of exterior radius 7?, interior 
radius r — 0.67?, and eccentricity 5 = 0.227?. Two openings 
of width W = 0.2227? are attached opposite to each other at 
an angle 9 £ [0, n/2] relative to the line connecting the circle 
centers, (b) Phase space, parameterized by the impact pa- 
rameter s and the transverse momentum component sin a of 
trajectories at the exterior circle. The lower panels show the 
phase space for open billiards with 9 = (c) and 9 — n/2 
(d), for which we only record trajectories that are injected 
through the openings (until they leave the billiard again). 



area larger than a Planck cell). 

A representative model for mixed regular and chaotic 
classical dynamics is the two-dimensional annular billiard 
p8[ , Fig. 0(a), which consists of the region between two 
circles with radii R, r, and eccentricity 5. Two leads 
(openings) of width W are attached opposite to each 
other at an angle 9 with respect to the axis through 
the two circle centers. The phase space can be param- 
eterized by the impact parameter s and the transverse 
component of the momentum sin a (with a the angle of 
incidence) of trajectories that are reflected at the exte- 
rior circle. The phase space of the closed annular billiard 
for r — 0.6R, 5 — 0.22R is shown in Fig. |l|(b), which 
displays two regular whispering-gallery (WG) regions, a 
large regular island, neighboring satellite islands, and a 
chaotic sea. Figs. |l|(c) and (d) show the phase space of 
the open annular billiard, which only includes trajecto- 
ries of particles that are injected into the billiard through 
leads attached at 9 — and 9 = ir/2, respectively (the 
width of the leads is W = 0.2227?). The large island is 
well coupled to one of the openings for 6 — 0, and is 
completely decoupled from both openings at 6 = O.bir. 
The chaotic sea and WG regions are well-coupled for ar- 
bitrary position of the openings, while the coupling of 
the satellite islands depends on 9 and W. In this way 
one can select regions in phase space by varying 9 and 
W. 

We computed the dimensionless conductance T = 
tiPt and the shot-noise suppression factor T = 
(2/7V) trvt(l — t't) from the transmission matrix t (the 
dimension of this matrix is given by the number of chan- 
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FIG. 2: (a) Dimensionless conductance T(E) and shot-noise 
suppression factor J-(E) for leads attached at 9 — 0.5n. (b) 
Energy-averaged suppression factor as a function of lead po- 
sition 9. 



nels N in each lead) jy. The transmission matrix is 
obtained numerically by the method of recursive Green 
functions fl9|] , for which space is discretized on a square 
lattice. In terms of the lattice constant a, R = 144 a, 
r = 86.4 a, 6 = 31.7 a, and W = 32 a. Energy E will 
be measured in units H 2 /(2ma 2 ) and time in units of 
2ma 2 /h, with m the mass of the charge carriers. In 
these units the mean level spacing A ~ 0.0003. We 
will work in the energy window E £ (0.408,0.433), in 
which the Fermi wavelength Af — 9.5 a, resulting in 
N = Int(2W/A F ) = 6. 

The upper two panels of Fig. |^ show T and T, as a 
function of energy, for the case that the two leads are 
attached at 9 — ir/2. The dependence of the energy- 
averaged suppression factor T on 9 is shown in the lowest 
panel of Fig. ||. Our first observation is that for all values 
of 9, the suppression factor is smaller than the universal 
value 1/4 for fully chaotic motion. Moreover, the shot 
noise is more suppressed for the cases that one opening 
couples to the large regular island (0 < 9 < 0.37r) than for 
the cases in which the regular islands is decoupled from 
the openings (9 > 0.37r). This behavior indicates that 
electrons injected into the large regular island contribute 
less to the shot noise. 

In order to determine in detail which processes are re- 
sponsible for the additional shot-noise suppression, we 
introduce a function /(t av ) which can be interpreted as 
the probability distribution function of deterministic pro- 
cesses with dwell time £dweii — t av . This will be achieved 
by monitoring the rate of change of the shot noise, 



/(*av) = -(dMav^^avV-^ch, 



(1) 



as we successively replace the dynamics on time scales 
?i t av by fully chaotic dynamics (equivalent to diffractive 
impurity scattering) . The function / is normalized to the 
total shot-noise suppression J°° di av /(i av ) = 1 — T jT Q ^, 
taken relative to the universal value of chaotic dynamics. 
The transmission matrix t determining T and J- is & 
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subblock of the scattering matrix 

where r and r' are reflection coefficients and the trans- 
mission matrix t' contains the same information as t. The 
elimination of the system-specific details on time scales 
> i av is achieved by averaging the scattering matrix S 
over an energy range [E — E av /2, E + E av /2] of width 
E av — h/t av (which will be taken inside the total energy 
range [0.408,0.433] of our numerical simulation), 

S(E av ;E ) = E- v 1 dES(E). (3) 

JE -E av /2 

Here the information on processes with longer time scales 
than i av is lost, because it is encoded into the short- 
range energy correlations (fluctuations) of the scatter- 
ing matrix |2fj| , while the information on the dynamics 
on shorter time scales than i av modulates the scatter- 
ing matrix on larger energy scales and hence is retained. 
Chaotic processes are then introduced to substitute the 
eliminated ones by coupling to an auxiliary chaotic sys- 
tem with scattering matrix So, taken from the appropri- 
ate circular ensemble of random matrix theory (observing 
the same symmetries as the original scattering matrix, as 
time-reversal or spatial parities pl|]), resulting in 

S'(E av ; E ; S ) = S(E av ; E ) + T'(l - SqK^SoT. (4) 

The ensemble of scattering matrices (Q) is the so-called 
Poisson kernel ^, |l5[ [l7j , with S the so-called optical 
scattering matrix. The coupling matrices T, T', and 1Z 
must be chosen such that S' is a unitary matrix, but 
the invariance of the circular ensemble guarantees that 
results do not depend on their specific choice. 

We now can calculate the mean suppression factor 
F(E av ) for fixed E av (and hence fixed i av ), first by 
averaging the noise within each Poisson kernel (fixing 
also Eq), and then averaging these values over Eq 6 
(0.408,0.433). The result for different positions of the 
leads is shown in Fig. 0. The value J-(E av = 0) is iden- 
tical to the suppression factor T shown in Fig. ||(b), 
because the Poisson kernel becomes a delta function at 
S(Eq) as S approaches this unitary matrix for E av — > 0. 
For increasing E av the suppression factor T(E av ) ap- 
proaches the universal value J- C h of random-matrix the- 
ory, because S = in this limit. The function ^{E^v) 
is monotonically increasing almost everywhere (negative 
values of the slope are within the statistical uncertainty in 
the numerical simulation) , indicating that the noise is en- 
hanced by the replacing the system-specific deterministic 
details of the transport by the indeterministic transport 
of random matrix theory. Consequentially, as shown in 
Fig. ||(b), the function /(i av ) introduced in Eq. (|j) is pos- 
itive and hence allows the interpretation of a probability 
distribution function of deterministic scattering. 




FIG. 3: (a) Shot-noise suppression factor T as a function of 
the strength of diffractive impurity scattering with a rate ~ 
E av /h, for different positions 9 of the leads, (b) Distribution 
function /(i a v) of deterministic processes for 9 — and 9 = 
0.5tv. 

The most striking feature in Fig. |^ is that F(E av ) rises 
quickly in the three cases 9 = 0, 0.1 7r, 0.2 it, in which the 
leads couple to the large regular island in phase space, 
while the slope of J-(E av ) almost vanishes up to an en- 
ergy £indet ~ 0.002 for 9 = 0.4 it, 0.5 tt (where the regu- 
lar region is decoupled). The distribution function f(t av ) 
decays only slowly in the former case, while vanishes al- 
most identically beyond a time i; n det ~ 500 in the latter 
case. This difference can be interpreted as follows. Wave 
packets are trapped for long times in regular regions, but 
show only little dispersion because of the stability of regu- 
lar trajectories p2| . Hence the scattering is deterministic 
and strongly affected by introducing diffractive scattering 
even at these long time scales. On the other hand, wave 
packets disperse quickly in the chaotic regions, where the 
transport is already indeterministic for i av > tindct, and 
this is not modified by adding diffractive scattering. The 
shot noise suppression for 9 = 0.4 n and 9 = 0.5tt hence 
arises from short chaotic trajectories with dwell times 
^ £indet- This is in accordance with the predictions of 
Rcf. JTo| ; chaotic quantum transport is deterministic up 
to the Ehrcnfest time (1/A) ln(27ri?/Ap) ~ 500, estimated 
from the Lyapunov exponent A ~ «f / R and the Fermi ve- 
locity v-p- Thus /(iav) reveals system-specific time scales 
as the Ehrcnfest time in chaotic regions and long dwell 
times in regular regions. Moreover, /(t av ) vanishes for 
iav ^ rnin idweii — 50 (not resolved in Fig. ||) . 

Our interpretation is confirmed by a semiclassical es- 
timate which discriminates the contributions to the shot 
noise from different regions in phase space. We divide 
the phase space into Planck cells of area 2W/ {NR) and 
calculate the fluctuations in the semiclassical occupation 
numbers of these cells. The average occupation number 
fj = TjlJl + Tjnfji is obtained from the classical trans- 
mission probabilities Tjl^ from cell j to the left (L) and 
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FIG. 4: Semiclassical estimates of T and / of Fig. [], obtained 
by introducing indeterminism into the classical dynamics be- 
yond a given dwell time tdweli = ^/^dwell- 



the hierarchical structure of small regular regions embed- 
ded in the chaotic part of phase space. A pronounced 
shot noise suppression should arise once that these trap- 
ping regions are resolved, because in the vicinity of stable 
structures the Ehrenfest time depends only algebraically 
on h [§2j. 

In conclusion, we have demonstrated that shot noise 
is a measure of the amount of deterministic transport 
through generic quantum dots. The dynamical analysis 
developed in this work may be used to study time scales 
in a large variety of systems. Experimentally, these time 
scales could be probed by measuring shot noise while 
tuning the indctcrministic scattering rate by adjusting a 
gate voltage Q . 

We thank S. Tomsovic for useful discussions and P. W. 
Brouwer for advice on implementing the Poisson kernel. 



right (R) opening, respectively, where Jl,r is the Fermi 
distribution in the electronic reservoirs attached to each 
opening. Within the model of minimal correlations of 
Ref. JtJ , the shot-noise suppression factor is then approx- 
imated as 



=</,•(!-/,•)}, 



(5) 



where (• ■ •) denotes the average over the phase space cells 
at the openings. 

Diffractive scattering can be introduced into the semi- 
classical theory by replacing the classical transmission 
probability of trajectories longer than a given dwell time 
^dwcii by the value 1/2 of ergodic dynamics. [The uni- 
versal value of shot noise for chaotic dynamics follows 
from Eq. (||) by inserting the ergodic value f rj — 1/2.] 
Figure | shows the resulting T c \(E Awc \\) and / c i(£dwcii), 
with -Edweii = ^Adweii- The same relative time scales 
for the different positions of the leads are found as in 
the exact results in Fig. ||. The agreement is reasonable 
given that we are not deep in the semiclassical limit and 
accounting for the fact that the model of minimal cor- 
relations in Ref. |?j originally was devised for chaotic 
dynamics. The suppression factor can be further de- 
composed into contributions of different regions in phase 
space. For 9 = 0, .T^i (indwell) initially increases due to 
the large regular island and then due to deterministic 
processes in the chaotic sea, while for 6 = 0.5tt it is con- 
stant up to -Edweii — -Eindot- Unlike the large regular 
island, the regular WG regions (which we have ignored 
so far in our discussion) contribute to shot noise with the 
universal value, which is explained by the fact that for 
the present system parameters Planck cells do not yet 
resolve reflected from transmitted trajectories which are 
injected into these regions. One can see from the model 
of minimal correlations, Eq. (||), that shot noise even- 
tually will be suppressed also in the WG regions when 
the semiclassical limit is approached. The same holds for 
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